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ABSTRACT 

Several  experiments  showed  that  the  classical  Kolmogorov  power  spectral  density  of  the  refractive-index  sometimes 
does  not  properly  describe  the  statistics  of  the  atmosphere.  In  this  paper  we  show  an  experimental  testbed  able  to 
generate  non-classical  Kolmogorov  turbulence  by  using  a  liquid  crystal  spatial  light  modulator.  The  testbed  is  used  at 
Naval  Postgraduate  School  for  laboratory  investigation  of  laser  beam  propagation  in  maritime  environment  where  a 
power  law  different  from  classical  Kolmogorov,  11/3,  could  be  present.  Applications  of  this  testbed  are  ship  to-ship  free 
space  optical  communication,  imaging  and  high  energy  laser  weapons 
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1.  INTRODUCTION 


Atmospheric  turbulence  is  one  of  the  most  serious  drawbacks  to  laser  beam  propagation  in  free  space  [1].  As  a  laser 
beam  begins  to  propagate  through  the  Earth’s  atmosphere,  the  varying  index  of  refraction  caused  by  the  turbulent  air  will 
alter  the  optical  path  of  the  light.  The  Earth's  atmosphere  can  be  described  as  a  locally  homogeneous  medium  in  which 
its  properties  vary  with  respect  to  temperature,  pressure,  wind  velocities,  humidity  and  other  factors.  Also,  the  Earth's 
atmosphere  temporally  changes  in  a  quasi-random  fashion  and  therefore  it  acts  like  many  random  lenses  present  in  the 
air.  All  of  these  processes  are  usually  simply  referred  to  as  “atmospheric  turbulence.”  The  Kolmogorov  model  of  energy 
distribution  in  a  turbulent  medium  is  a  useful  statistical  model  to  describe  the  fluctuation  of  the  index  of  refraction  due  to 
temperature  changes  [1].  This  model  is  supported  by  a  variety  of  experimental  measurements  and  is  quite  widely  used  in 
simulations  for  the  propagation  of  electromagnetic  waves  through  a  random  medium.  Although  for  many  years, 
scientific  communities  have  investigated  optical  turbulence  present  in  the  atmosphere  and  new  important  results  have 
been  carried  out,  the  atmosphere’s  statistical  behavior  is  still  not  well  understood.  For  example,  the  power  spectrum  of 
the  index  of  refraction  sometimes  does  not  follow  the  classical  Kolmogorov  power  law,  11/3,  and  it  rather  follows 
different  power  laws  [2-4].  NRL  researchers  recently  measured  scintillation  in  maritime  environment  [9],  some 
conclusion  are  that  classical  Rytov  approximation  [1]  might  not  be  representative  of  propagation  in  a  maritime 
environment,  scintillation  shows  a  much  faster  transition  into  strong  saturation  regime,  refractive  index  structure 
parameter,  which  is  a  measure  of  the  strength  of  the  fluctuations  in  the  refractive  index  in  the  atmosphere,  is  more  stable 
during  the  day,  and,  more  importantly,  refractive  index  power  spectra  is  sometimes  different  from  Kolmogorov  [11- 
12]. The  Navy  has  an  interest  in  a  laser-based  directed  energy  weapon  (DEW)  system  because  of  the  need  for  anti-ship 
cruise  missile  and  tactical  air  defense  [5-6].  US  Naval  Research  Laboratory  has  undertaken  a  multi-year  research  effort 
in  free  space  optical  communications  in  maritime  environment  [7-12].  For  these  and  other  reasons  is  important  and 
strategic  to  do  more  investigation  on  laser  propagation  through  atmospheric  turbulence  in  a  maritime  environment  [17], 
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where  Kolmogorov  power  spectrum  does  not  describe  well  the  statistical  behavior  of  the  index  of  refraction  fluctuations 
in  the  atmosphere. 

The  refractive  index  structure  parameter  is  known  to  behave  differently  in  marine  and  terrestrial  environments,  since  a 
large  body  of  water  acts  quite  differently  from  soils  [13].  Friehe  [14]  investigated  the  influence  of  humidity  on  the 
optical  refractive  index  in  the  maritime  environment,  and  evidence  of  non-Kolmogorov  power  spectrum  is  even  shown. 
Finally,  significant  evidence  of  power  law  deviation  from  classical  Kolmogorov  turbulence  is  shown  in  [15],  where  is 
investigated  the  effect  of  intermittency  in  fully  developed  turbulence  on  the  spectra  of  fluctuations  of  “passive  scalars”. 

In  this  paper  we  describe  a  testbed  used  at  Naval  Postgraduate  School  for  laboratory  investigation  of  laser  beam 
propagation  in  maritime  environment  where  a  power  law  different  from  classical  Kolmogorov,  11/3,  could  be  present. 
Applications  of  this  testbed  are,  for  instance,  ship  to-ship  free  space  optical  communication,  imaging  and  high  energy 
laser  weapons. 


2.  NON-CLASSICAL  KOLOMOGOROV  POWER  SPECTRUM 


The  traditional  Kolmogorov  spectrum  is  primarily  used  in  theoretical  studies  due  to  its  simplicity  and  it  is  only  valid  in 
the  inertial  subrange,  where  no  inner  scale  and  outer  scale  effects  are  included. 


In  this  paper  we  use  the  power  spectrum  introduced  in  [4]  with  a  power  law  OC  that  can  deviate  from  the  classical 
Kolmogorov.  The  non-Kolmogorov  power  spectrum  is  given  by 


On(K,a)  =  A(a)-Cl- 


0<k<oo,  2<a<4 


(l) 


where  a,  is  the  spectral  index  or  power  law,  C  ;=p<  is  a  generalized  structure  parameter  with  units  m 
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KQ  = - ,  is  the  outer  scale  parameter; 

L0 


K  =^Ui 


is  the  inner  scale  parameter  ; 


and  the  symbol  T  j  is  the  gamma  function. 
The  constant  is  defined  by 

c(a)  =  {|K./l(a)-r[^j}a5 


(2) 


(3) 


The  procedures  to  carry  out  and  C^Ot^  are  shown  in  Andrews  and  Phillips  [1],  pag. 74-75,  or  in  [4]  where 

essentially  the  asymptotic  behavior  of  the  structure  function  has  been  considered. 
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To  show  the  non-Kolmogorov  power  law,  we  plot  in  Fig.l  the  power  spectrum  as  a  function  of  wavenumber  K  for  three 
different  values  of  CL ,  we  set  l0  =  0.0\m-,L0  =  \0m,C2n=  1.7-10  "14m3“a. 


Fig.  1  -Power  spectrum  as  a  function  of  wavenumber  for  three  different  power  law  values 

We  deduce  from  Fig.  1  that  a  different  power  law  from  the  classic  value  of  Kolmogorov,  CL  =  1 1  /  3,  changes  the  slope 
of  the  spectrum  and  obviously  impacts  on  atmospheric  turbulence  parameters  such  as  phase  structure  function, 
scintillation  and  long  term  beam  spread  [17].  Therefore  applications  of  laser  beam  propagation  in  maritime  environment, 
such  as  ship-to-ship  free  space  optical  communication  and  imaging  should  consider  this  non-Kolmogorov  effect  into  the 
budget  link  analysis. 


3.  TURBULENCE  GENERATION  METHOD 


The  Naval  Research  Laboratory  has  developed  a  new  method  for  generating  atmospheric  turbulence  and  a  testbed  that 
simulates  its  aberrations  far  more  inexpensively  and  with  greater  fidelity  than  many  other  methods  [16].  This  system 
uses  a  Liquid  Crystal  (LC)  Spatial  Light  Modulator  (SLM)  which  allows  the  simulation  of  atmospheric  seeing  conditions 
ranging  from  very  poor  to  very  good  and  different  algorithms  may  be  easily  employed  on  the  device  for  comparison. 
These  simulations  can  be  dynamically  generated  and  modified  very  quickly,  easily  and  multiple  devices  can  be  used  to 
simulate  multiple  layers  of  turbulence  in  a  laboratory  environment.  Current  testing  with  using  multiple  LC  SLMs  is 
under  investigation  at  the  Naval  Research  Laboratory  and  the  Naval  Postgraduate  School. 

The  method  described  in  [16]  permits  to  generate  the  classical  Kolmogorov  atmospheric  turbulence.  In  this  paper  we 
extend  the  method  even  for  non-Kolmogorov  turbulence.  However  here  we  use  Zemike  polynomials  instead  than 
Karhunen-Loeve  polynomials,  therefore  the  coefficients  will  not  be  statistical  independent.  A  future  step  will  be  to  carry 
out  and  implement  the  coefficients  for  phase  expansion  by  using  Karhunen-Loeve  polynomials. 

It  is  well  known  that  Zernike  polynomials  represent  aberrations  from  low  to  high  order  with  increasing  order  index.  A 
wavefront  can  generally  be  represented  by 
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(4) 
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where  the  CV.  s  are  the  amplitudes  of  the  aberrations  for  the  jth  Zernike  order,  Z . ,  and  M  is  the  total  number  of  orders 
in  which  the  wavefront  is  represented. 

Noll  [18]  analyzed  the  phase  distortions  due  to  the  atmosphere  in  terms  of  Zernike  polynomials  for  Kolmogorov 

turbulence;  Boreman  and  Dainty  [19]  described  how  to  calculate  in  closed  form  the  CL\ s  coefficients  of  the  Zernike 

expansions  for  non-Kolmogorov  turbulence  using  the  power  spectrum  (1).  As  main  result,  for  zero  inner  scale  and 
infinite  outer  scale,  they  carry  out  that 
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valid  for  n  >  1  and  3  <  a  <  4.  In  expression  (5)  il  is  the  radial  order  as  shown  in  (16),  D  is  the  aperture  diameter, 
and  is  the  Fried  parameter  as  defined  in  (17). 


Also,  Boreman  and  Dainty  show  in  [19]  that  the  CL\  s  coefficients  of  the  Zernike  expansions  for  non-Kolmogorov 

turbulence  using  power  spectrum  (1)  and  including  the  effect  of  inner  and  outer  scale  can  be  numerically  obtained  by 
numerically  solving 
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However,  in  this  paper  we  will  use  only  equation  (5),  therefore  we  will  not  consider  inner  and  outer  scale  effects. 

It  is  well  known  that  the  temporal  phase  variances  of  atmospheric  turbulence  have  a  Gaussian  random  distribution, 
therefore,  by  taking  Equation  (4)  and  modifying  it  such  that  there  is  Gaussian  random  noise  factored  in  gives, 

M 

®  =  I 

7=1 

where,  X .  is  the  amount  of  noise  for  the  jth  mode  based  on  a  zero-mean  unitary  Gaussian  random  distribution  and  it 
can  be  generated  by  using  a  software  random  number  generator  that  uses  the  Gaussian  distribution.  However  generating 
a  continuous  transition  for  the  atmospheric  turbulence  realization  temporally  will  require  another  method.  The  X’.s  can 
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be  modified  from  being  just  random  numbers  to  a  continuous  function  of  time  for  each  mode.  Thus,  Equation  (7)  can  be 
rewritten  as 


unitary  Gaussian  distribution.  Next,  a  spline  curve  is  fit  to  this  vector,  as  shown  in  Fig.  2,  and  this  spline  curve  is  now 

the  X  .{t^  temporal  function  for  generating  the  wavefronts  in  the  atmospheric  turbulence  simulation.  This  process  is 

repeated  for  each  of  the  M  modes  in  the  model.  Without  this  splining  technique,  the  change  between  phase  screens 
would  be  discontinuous  and  would  not  provide  an  accurate  representation  of  the  atmosphere  for  testing  an  adaptive 
optics  system.  In  reality,  the  Earth’s  atmosphere  is  a  continuous  medium.  With  this  Spline  technique,  the  temporal 
transition  of  the  wavefronts  in  the  atmospheric  turbulence  simulation  is  continuous  and  smooth. 


Fig.2-Sample  Xi(t)  temporal  function  generated  from  a  vector  of  a  few  Gaussian  random  number 


Finally,  Wilcox  [16]  showed  that  the  method  of  generating  phase  distortions  for  atmospheric  turbulence  simulation  via  a 
Spline  technique  is  virtually  the  same  as  the  Frozen  Seeing  method  with  the  advantage  of  being  far  less  computationally 
intensive  on  a  computer  system.  A  comparison  between  the  new  Spline  method  and  the  Frozen  Seeing  method  for 
temporal  simulation  of  atmospheric  turbulence  was  performed  [16].  After  generating  and  analyzing  10  different 
realizations  of  atmosphere  from  the  Frozen  Seeing  method  and  performing  the  cross-correlation  between  the  two 
methods’  temporal  transitions  it  was  found  that  the  two  methods  are  similar  and  consistent  with  each  other. 

4.  SLM  CONTROL  SOFTWARE 


While  the  diffractive  nature  of  the  SFM  decreases  the  available  light  and  provides  additional  alignment  challenges,  the 
device  provides  a  great  deal  of  flexibility  in  generating  and  applying  aberrations  in  the  optical  path.  Atmospheric 
scenarios  can  be  changed  quickly  in  software  without  having  to  change  hardware.  This  provides  a  significant  advantage 
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over  other  popular  hardware  such  as  rotating  plates  imprinted  with  specific  atmospheric  statistics,  which  cannot  be 
changed  without  constructing  new  plates. 

The  Naval  Research  Laboratory  (NRL)  has  developed  software  to  apply  atmospheric  aberrations  on  the  SLMs  using  a 
Matlab  graphical  user  interface  (GUI).  The  SLM  control  GUI  is  shown  in  Fig.3.  The  atmosphere  generated  in  software 
is  based  on  non-Kolmogorov  statistics  as  previously  described.  As  analytical  models  are  developed  to  describe 
horizontal  turbulence  and  include  effects  particular  to  the  maritime  environment,  the  software  can  be  modified  to 
accommodate  these  model  changes. 

The  NRL  software  GUI  allows  the  user  to  specify  telescope  and  site  parameters  such  as  the  telescope  diameter,  D ,  and 
the  atmospheric  coherence  length  or  Fried  parameter,  r0  .  Non-Kolmogorov  turbulence  generation  can  be  performed 

with  the  Zemike  polynomial  expansion  (7),  the  Karhunen-Loeve  expansion  can  be  used  only  for  Kolmogorov 
turbulence.  The  Karhunen-Loeve  modes  contain  a  statistically  independent  set  of  coefficients  based  on  Zernike  modes, 
and  are  often  used  in  turbulence  simulation.  Next  step  will  be  to  extend  the  non-Kolmogorov  turbulence  generation  even 
for  the  Karhunen-Loeve  expansion. 

Once  site  parameters  and  any  initial  alignment  biases  or  static  aberrations  are  selected,  the  controller  generates  a  user- 
specified  number  of  phase  screens  to  represent  the  atmospheric  simulation.  At  the  present,  the  value  of  non-Kolmogorov 
power  law,  OC5  is  not  shown  on  the  GUI,  but  it  is  set  inside  a  Matlab  function  which  is  part  of  the  code.  In  Fig.  4  an 
example  of  phase  screen  is  shown.  A  secondary  annular  obscuration  can  be  included  in  the  simulation  to  simulate  a 
telescope’s  secondary  mirror.  By  pushing  the  button  “Generate  Turbulence”  on  the  GUI  several  frames  of  turbulence 
are  generated  as  shown  in  the  Video  reachable  at  http://dx.doi.Org/10.l  1 17/12.903681.1. 

The  simulation  is  then  applied  on  the  ferroelectric  Spatial  Light  Modulator,  from  Boulder  Nonlinear  Systems,  which  is 
able  to  operate  up  to  1  kHz,  it  has  512x512  pixels  and  an  active  area  of  7.7mm  x  7.7mm. 

The  simulation  can  be  changed  easily  and  quickly  as  new  atmospheric  and  site  parameters  are  desired  for  testing,  and  the 
underlying  code  can  be  changed  as  new  models  of  turbulence  become  available. 


Fig.  3  -  Matlab  graphical  user  interface  developed  by  NRL 
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Fig.  4  -  An  example  of  phase  screen  applied  on  BNS  SLM 

5.  EXPERIMENTAL  TESTBED 


An  experimental  setup,  as  shown  in  Fig.  5,  has  been  built  at  Naval  Postgraduate  School  to  conduct  the  practical  tests. 


Fig.  5  -  Experimental  setup  built  at  Naval  Postgraduate  School 
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The  system  is  configured  with  two  parallel  beams  that,  after  coming  through  a  microscope  from  the  laser,  are  split  by  a 
polarized  beam  splitter,  PBS1.  One  of  these  beams  is  the  reference  beam  that  is  propagated  until  the  camera  by  using  a 
mirror  and  a  beam  splitter,  BS3.  The  Spatial  Light  Modulator  is  a  ferroelectric  from  Boulder  Nonlinear  Systems  with 
512  x  512  pixels  and  an  active  area  of  7.7mm  x  7.7mm.  The  second  beam  going  from  PBS1  goes  through  another  beam 
splitter  PBS2  and  is  directed  onto  the  SLM.  After  being  reflected  from  the  SLM  the  beam  comes  back  through  PBS2  and 
into  a  couple  of  lenses  and  an  iris  placed  at  the  focus  to  select  only  the  first  diffraction  order  of  the  SLM.  It  then  travels 
to  BS4,  where  it  interferes  with  the  reference  beam.  A  tip-tilt  bias  on  the  phase  was  introduced  to  “send”  light  to  the  first 
diffraction  order.  The  amplitude  and  phase  information  is  carried  by  the  first  diffraction  order  selected  by  the  iris. 

6.  CONCLUSION 


In  this  paper  we  described  a  method  of  using  a  Spatial  Light  Modulator  to  simulate  non-Kolmogorov  turbulence  at 
laboratory.  The  method  of  generating  phase  distortions  for  atmospheric  turbulence  simulation  via  a  Spline  technique  is 
virtually  the  same  as  the  Frozen  Seeing  method  with  the  advantage  of  being  far  less  computationally  intensive  on  a 
computer  system.  This  advantage  can  be  exploited  in  the  development  of  a  software  package  that  can  drive  any  LC  SLM 
to  simulate  atmospheric  turbulence  at  almost  any  wavelength  for  any  telescope  diameter  at  any  location.  AO  and  Free 
Space  Laser  Communication  systems  can  also  be  tested  for  performance  evaluations. 
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